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ABSTRACT 



We investigate the soft X-ray transients with black hole primaries which 
may have been the sources of gamma-ray bursts and hypernovae earlier in 
their evolution. For systems with evolved donors we are able to reconstruct 
the pre-explosion periods and find that the black-hole mass increases with the 
orbital period of the binary. This correlation can be understood in terms of 
angular-momentum support in the helium star progenitor of the black hole, if 
the systems with shorter periods had more rapidly rotating primaries prior to 
their explosion; centrifugal support will then prevent more of its mass from 
collapsing into the black hole on a dynamical time. This trend of more rapidly 
rotating stars in closer binaries is usual in close binaries, and in the present 
case can be understood in terms of spin-up during spiral in and subsequent 
tidal coupling. We investigate the relation quantitatively and obtain reasonable 
agreement with the observed mass-period correlation. An important ingredient 
is the fact that the rapidly rotating new black hole powers both a GRB and the 
hypernova explosion of the remaining envelope, so that the material initially 
prevented from falling into the black hole will be expelled rather than accreted. 

For systems in which the donor is now, and will remain, in main sequence 
we cannot reconstruct the pre-explosion period in detail, because some of their 
history has been erased by angular momentum loss through magnetic braking 
and gravitational waves. We can, however, show that their periods at the 
time of black hole formation were most likely 0.4-0.7 day, somewhat greater 
than their present periods. Furthermore, their black holes would have been 
expected to accrete ~ IM© of material from the donor during their previous 
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evolution. Comparison with predictions suggests that httle mass will be lost in 
the explosion for the relatively high pre-explosion periods of these binaries. 

A natural consequence of the He star rotation is that black holes formed in 
the shorter period (before explosion) soft X-ray transients acquire significant 
Kerr parameters. This makes them good sources of power for gamma-ray bursts 
and hypernovae, via the Blandford-Znajek mechanism, and thus supports our 
model for the origin of gamma-ray bursts in soft X-ray transients. 

Subject headings: black hole physics — stars: binaries: close — accretion — 
gamma-ray bursts 
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1. Introduction 

Recent observations strongly suggest a connection between gamma-ray bursts and 
supernovae, with indication that the supernovae in question are especially energetic and 
of type Ib/c, i.e., core collapses of massive stars which have lost their hydrogen envelope 
(see Van Paradijs, Kouvehotou, and Wijers 2000, and references therein). This supports 
suggestions by Woosley (1993) and Paczynski (1998) for the origin of gamma-ray bursts 
in stellar core collapses. The hydrodynamics of a jet escaping from a star and causing 
its explosion was explored in detail by MacFadyen and Woosley (1999), who showed that 
contrary to accepted wisdom, a fairly baryon-free, ultra-relativistic jet could plow through 
the collapsing star and emerge with large Lorentz factors. The powering of the outflow 
by coupling of high magnetic fields to the rotation of the black hole (Blandford & Znajek 
1977), first suggested by Paczynski (1998) in the context of gamma-ray bursts, was worked 
out in detail by Van Putten (1999, 2001). Li has also discussed the deposition of energy 
from a black hole into the accretion disk in a recent series of papers (e.g., Li 2000a). 

Building on these thoughts, we have modeled both the powering of a gamma-ray burst 
by black-hole rotation and the stellar evolution pathways that set up favorable conditions 
for that mechanism (Brown et al. 2000). An essential ingredient in this model is a rapidly 
rotating black hole, and it is this aspect that we focus on in the present paper. A single 
star initially in uniform rotation will tend to develop a differential rotation because the core 
contracts strongly during evolution, and angular-momentum conservation will therefore 
increase its angular velocity. However, given enough time viscous stresses will even out 
these differences, and thus the net result is a loss of angular momentum of the innermost 
regions of the star. Spruit & Phinney (1998) argued that magnetic-field mediated coupling 
is strong enough in single stars during the giant phase to make the cores very slow; so 
slow, in fact, that they required an asymmetric kick in the birth of pulsars to get their spin 
frequencies up to observed values. Livio & Pringle (1998) subsequently used observations 
of novae to argue for a weaker coupling, but still their coupling strength would lead to spin 
energies of black holes that are negligible as power sources for gamma-ray bursts. 

However, as suggested by MacFadyen & Woosley (1999), a massive star in a close 
binary will spin faster for a number of reasons: first, when the hydrogen envelope is lifted 
off by spiral-in, it will cease to serve as a sink of angular momentum for the core. Second, 
the tidal friction concomitant of the spiral-in process will spin up the inner region, giving it 
a larger angular momentum than the same region in a single star (Rasio and Livio 1996). 
Third, tidal coupling in the close binary will tend to bring the primary into corotation with 
the orbital period. This latter process is not very efficient in the short post spiral-in life 
of the binaries we consider, but its effect does probably matter to the outer layers of the 
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helium star, which can be important for our work. With its more rapid rotation, the hehum 
star then forms a black hole with a large Kerr parameter, which immediately after its 
formation (in a few seconds) begins to input power into its surroundings at a very high rate. 
This, then, powers both a gamma-ray burst (e.g.. Brown et al. 2000) and the expulsion of 
the material that was centrifugally prevented from falling into the black hole. In fact. Van 
Putten (1999, 2001) estimates that the power input into that material exceeds that into the 
GRB and Li (2000b) also finds that more energy can be extracted by the disk than by the 
GRB. It should be noted that an initially less rapidly rotating black hole could be spun up 
by disk accretion quite rapidly, and start a similar process after some accretion has taken 
place (MacFadyen & Woosley 1999, Brown et al. 2000). Some implications of such more 
complicated sequences of events are discussed by Lee, Lee, & Van Putten (2001). 

In section ^ we present the data on known soft X-ray transients, showing the relation 
between present orbital period and black-hole mass. Since theory predicts a relation 
between pre-explosion orbital period and black-hole mass, we then consider carefully the 
pre- and post-explosion evolution of the systems (section and use this to reconstruct the 
pre-explosion orbit for as many systems as possible (section ^). Then we develop our model 
for the mass and spin of black holes in soft X-ray Transients (SXTs), and use it to explain 
the mass-period correlation (section]^). We summarize our conclusions in section ^. 

2. An empirical correlation between orbital period and black-hole mass 

We have collected data from the literature on black-hole binaries in our Galaxy. In 
table 1 we collect data of those in which the mass function is known, and some manner of 
mass estimate for both the black hole and companion can be given. In table 2 we list the 
properties of two key systems in more detail. In Fig. |I| we show the masses of the black holes 
as a function of orbital period. While the ranges of black-hole masses for main-sequence 
and evolved systems overlap, the latter tend to have higher masses; the exception is Nova 
Scorpii 1994, which we shall see later is a natural but rare case of the general evolution 
scenario that we describe in this paper. In Fig. we show the donor masses as a function 
of orbital period. They show a more obvious trend of more massive donors in evolved 
systems. As we shall see, this is a natural consequence of the fact that only evolved systems 
can come into Roche contact in wide binaries, and more massive donors are more likely to 
come into contact via nuclear evolution. (The various curves are explained in section |^.) 

In the following sections, we shall argue that the correlation between black-hole mass 
and period also has physical meaning: the shorter the orbital period, the more rapidly 
rotating the helium star progenitor to the black hole. Rapid rotation centrifugally prevents 
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some fraction of the helium star to collapse into a black hole, resulting in a smaller black 
hole mass. The reason why the correlation in Fig. |1] is weak is that evolution of the 
binary since formation of the black hole has washed out the relation. Properly, we should 
consider the correlation between pre-explosion orbital period and post-explosion black hole 
mass. Much of our work presented here is concerned with understanding the evolution of 
these binaries, and using this knowledge to find the systems for which we can reconstruct 
those parameters. Using that subset, we show much better agreement between our model 
predictions and the observed relation between reconstructed period and mass; this supports 
our evolutionary model and has ramifications for the origin of gamma-ray bursts. 

3. The evolution of soft X-ray transients 

3.1. Prior to the formation of the black hole 

Following on the work of Brown, Weingartner, & Wijers (1996), who showed the 
importance of mass loss of helium stars in binaries in determining the final outcome of 
binary evolution. Brown, Lee, & Bethe (1999), Wellstein & Langer (1999), and Brown et 
al. (2001a) showed that massive helium stars could evolve into high-mass black holes only 
if they were covered with hydrogen during most of their helium core burning era (Case C 
mass transfer in binaries). In case A or B mass transfer in binaries (Roche Lobe overflow in 
main-sequence or red-giant stage) the Fe core that was left was too low in mass to go into 
a high-mass black hole. Brown et al. (2001a) showed that high-mass black holes could be 
formed only if the mass was taken off of the black hole progenitor after helium core burning 
was finished; i.e. Case C mass transfer. Brown, Lee, & Tauris (2001c) showed that with 
the Schaller et al. evolution, this could happen only in the neighborhood of ZAMS mass 
2OM0, definitely not at 25Mq and higher: Because of the wind losses, the mass transfer 
would begin as Roche lobe overflow only in Case B mass transfer, with these higher main 
sequence masses, because the giant radii of these stars exceed their radii in the supergiant 
phase. Below we shall see that the high black- hole masses in a few binaries require us to 
extend the range upward to about 30 Mq. 

A major uncertainty in the evolution of all compact X-ray binaries is the phase of 
spiral-in that occurred in their evolution: these binaries are initially very wide, and when 
the primary fills its Roche lobe and transfers mass to the secondary, the mass transfer 
leads to instability, resulting in the secondary plunging into the primary's envelope. Next, 
dissipation of orbital energy of the secondary causes the primary's envelope to be ejected, 
and the orbit to shrink. Following the original work by Webbink (1984), Brown, Lee, & 
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Tauris (2001c) write the standard formula for common envelope evolution as 

XR XvLtti \ 2a/ 2ai J ^ ' 

where Mp is the total mass of the BH progenitor star just before the common envelope 
forms, Me is the mass of its hydrogen envelope, Mne is the mass of its core, Oi and Of are 
the initial and final separation, before and after the common envelope, respectively, and 
tl = Rl/cl is the dimensionless Roche- lobe radius. This equation essentially relates the 
loss of orbital energy of the secondary to the binding energy of the ejected envelope. The 
parameter A is a shape parameter for the density profile of the envelope. It can vary greatly 
between stars (Tauris & Dewi 2001), but for the extended, deeply convective giants we deal 
with in case C mass transfer it is always close to 7/6. (See also Appendix C of Brown et 
al. 2001b.) The parameter ctce accounts for the efficiency with which orbital energy is used 
to expel the envelope, and may also account for some other effects such as extra energy 
sources and the possibility that each mass element of the envelope receives more than the 
minimum energy needed to escape (see, e.g., Bhattacharya & Van den Heuvel 1991 and 
references therein). 

Given the parameters of the system at first Roche contact, when spiral-in starts, the 
final separation is determined by the product of A and ctce, the efficiency of the energy 
conversion. In general, these parameters are only the simplest recipe prescription for the 
complex hydrodynamical interaction during spiral-in. While we therefore cannot predict 
the value of Accce from first principles, we can try to find its value from constraints in some 
systems, and then assume it is the same for all similar systems. Brown, Lee, & Tauris found 
a great regularity in the evolution of SXTs with main-sequence companions, all but one of 
which are K or M stars, which constrained the efficiency Xa^e to be 0.2-0.5. However, these 
authors did not include mass loss in the explosion, which we shall do here in our evolution of 
SXTs with evolved companions. Since mass loss substantially widens the orbits, including 
it the common envelope evolution must bring the (pre-explosion) a/ to a smaller value: if 
Mpost is the black hole plus companion mass, and AM the mass lost in the formation of the 
black hole, we have 

,pre {1 + AM /Mpost)- (2) 

(after the orbit has been re-circularized). Therefore, there has been an extra widening since 
the explosion by a factor of up to about 1.5. We found Actce fo be in the lower part of the 
interval found by Brown, Lee, & Tauris (2001c), Xace ~ 0.2. 

We can achieve a more precise 'calibration' of the value of Actce if we manage to find 
some systems in which we can estimate both the initial and final separation. To estimate 
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the initial separation (at the onset of spiral-in) we need to know the mass and radius of 
the black-hole progenitor and combine this with the Roche-lobe filling condition. The 
helium star progenitors in at least three of the evolved binaries seem to be too massive 
for the 20 — 23Mq ZAMS progenitors used by Brown, Lee, & Tauris (2001c): the black 
hole in V404 Cyg is probably at least IOMq (Shahbaz et al. 1994, Shahbaz et al. 1996, 
Bailyn et al. 1998), and the black hole in Nova Scorpii is of mass ~ 5.4 ± O.SMq (Beer & 
Podsiadlowski 2001) and the mass loss in black hole formation is ^SMq (Nelemans, Tauris, 
Van den Heuvel 1999) so that the progenitor of the helium star must have been ~ IIMq. 
From Table |, the black hole in V4641 Sgr is of mass 9.6llo:88M0 (Orosz et al. 2001). The 
tentative conclusion from the above is that at least these binaries with evolved companions 
seem to have come from helium cores of ~ HMq, or ZAMS mass ~ 30Mq. With high wind 
mass loss rates as proposed by Schaller et al. (1992), such massive stars have larger radii as 
giants than as supergiants, thus making case C mass transfer impossible. However, since 
radii and mass loss rates of evolved stars are very uncertain, we take the view that the 
need for ~ IIMq helium cores implies that their progenitors, 30Mq main-sequence stars, 
do expand enough to allow case C mass transfer. 

In Fig. ^, we summarize the radius at the end of the giant branch as a function of 
ZAMS mass (Schaller et al. 1992). The ZAMS mass dependence of this final giant radius is 
adequately represented by a linear function in the region of 20 — 4OM0. We assume that 
the radial expansion during the helium burning can be scaled to the case of a 2OM0 star 
using this linear relation as follows: 

Further, we took the mass loss rate of 2OM0 as standard, and scaled the mass loss rate in 
proportion to the ZAMS mass. The allowed range of Case C mass transfer with ZAMS 
mass 2OM0 is 971i?0 < R < 1185i?0, that of Schaller et al. (1992). In Fig. | are given 
the possible initial orbital separations for Case C mass transfer for the I.9IM0 companion 
appropriate for Nova Scorpii (see section |^) and for the 6.53M0 companion appropriate for 
V4641Sgr (Orosz et al. 2001). 

Now if we look at eq. (|1|), we see that af scales almost linearly with the donor 
(companion) mass M^. The envelope mass Me is roughly O.TMgiant (Bethe & Brown 1998) 
and we use 

Mne = O.O8(M,,,„i/M0)i-^5M0 (4) 

so that 
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assuming \(Xce to be constant and with neglect of the small term in a~ in the r.h.s of 
eq. (0). From our curves, Fig. ^, we see that the 20% possible variation in results in 
the same percentage variation in a^. Because the actual ZAMS mass can be anywhere 
in the range 20 — SOM© there can be an additional ~ 25% variation in a/ with giant 
mass, as compared with the linear dependence on M^. In view of the modest size of these 
variations at a given donor mass, we make the approximation in the rest of the paper that 
the pre-explosion orbital separation depends only on M^, and scales linearly with M^. This 
simple scaling and the modest amount of scatter around it are partly the result of the weak 
dependences on initial parameters in eq. (^, but chiefly the result of the fact that our 
model uses case C mass transfer. This constrains the Roche contact to first occur when 
the radius of the star is in a very narrow range, between the maximum radius in the giant 
phase and the maximum radius in the supergiant phase. 

To complete our calibration of the spiral-in efficiency, we need to find systems in which 
we can also estimate the orbital separation just after spiral-in well. This is complicated by 
the fact that mass transfer has taken place since the spiral-in. Most SXTs have small mass 
ratios, and for such small mass ratios the orbital separation is fairly sensitive to the amount 
of mass transferred, making it hard to derive the post-spiral-in separation from the present 
one. The exception is V4641 Sgr, in which the present mass ratio is close to 1. Since the 
initial mass ratio could not have been significantly greater than 1 (since that would result 
in unstable mass transfer), and furthermore the orbital period changes very little with mass 
transfer for nearly equal masses, we can fairly approximate the post-spiral-in separation by 
the present one. In Fig. ^ we show the predicted ranges of post-spiral-in orbital periods for 
different values of \ace- Clearly, a value quite close to 0.2 is indicated. For 4U 1543— 47 (IL 
Lupi), we find that it is near the boundary between evolved and main-sequence evolution. 
To place it there, as discussed in section ^, we find from the reconstructed orbital period in 
Fig. ^ that Actce ~ 0.2 is also consistent with the properties of this system. 

In short, the general properties of SXTs and the specific cases of V4641 Sgr and IL 
Lupi favor Aoce ~ 0.2, which we therefore adopt as a general efficiency for the evolution of 
other transient sources. This then makes it possible to make quite specific predictions for 
the prior evolution of many of the other SXTs. 



From the above theory, certain regularities follow for the system behavior as a function 



3.2. Expected regularities 
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previous section). Furthermore, the relation between Roche lobe radius and donor mass 

1 /3 

when Md <C Mbh implies that RL/df oc M^' (e.g., Eggleton 1983). As a result, the Roche 
lobe radius of the donor just after spiral-in will scale with donor mass as Rl oc M^l On 
the other hand, the donor radius itself depends on its mass only as Rd oc M^'®. Therefore, a 
low-mass donor overfills its Roche lobe immediately after spiral-in. In the donor mass range 
we consider (M^ ^ 0.7Mq) it does not overfill its Roche lobe by much, so we assume that 
the system adjusts itself quickly by transfer of a small amount of mass to the He star, which 
widens the orbit until the donor fills its Roche lobe exactly. Above this minimum mass, 
there will be a range of donor masses that are close enough to filling their Roche lobes after 
spiral-in that they will be tidally locked and will come into contact via angular momentum 
loss (AML). Above this, there will be a range of mixed evolution, where both AML and 
nuclear evolution (Nu) play a role. Finally, for the most massive donors, > 2Mq, the 
post-spiral-in orbits will be too wide for AML to shrink them much, so mass transfer will 
be initiated only via nuclear expansion of the donor. Of course, the ranges of case C radii 
of stars and variations of primary masses will ensure that the boundaries between these 
regions are not sharp: near the boundaries the fate of the system depends on its precise 
initial parameters. 



4. Reconstructing the pre-explosion orbits 

Nova Sco 94 (GEO J1655—40): The most extensive evolutionary studies have been 
made for Nova Scorpii. Starting from the work of Regos, Tout & Wickramasinghe (1998) 
who make the case that the companion is in late main sequence evolution. Beer & 
Podsiadlowski (2001) carry out extensive numerical calculations of the evolution, starting 
with a pre-explosion mass of 2.5Mq and separation of ~ 6Rq. More schematically we 
arrived at a pre-explosion mass of I.9IM0 and separation of 5.33i?0. We consequently have 
an 0.4 day pre-explosion period. With ~ 6Mq mass loss in the explosion (Nelemans et al. 
1999), nearly half the system mass, the binary period increases to 1.5 day, well beyond the 
period gap. This is also the period required if the common-envelope efficiency in this binary 
was again 0.2 (Fig. This explains why Nova Sco is the only system with a giant donor 
and a black-hole mass in the lower end of the range: its evolution really places it among 
the narrow-orbit systems. Generally, the mass loss during explosion is mild, and does not 
change which category a system belongs to. But in those exceptional cases where the mass 
loss comes close to half the total mass, the orbit widens very much and converts an AML 
system to a nuclear-evolution system. We shall discuss in section |^ that help in expelling 
the mass may come from early onset of the GRB mechanism. After explosion the binary 
evolves to its present period by nearly conservative mass transfer. Our estimate is that 
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O.41M0 is transferred from the donor to the black hole. Brown et al. (1999) first made the 
case that Nova Scorpii was the relic of a GRB. 

V464I Sgr: As we discussed in section this system is our calibrator for the spiral- in 
efficiency, and we assume that its present state is very close to the one immediately following 
spiral-in. 

GRS 1915+105: Recently Greiner et al. (2001) have determined the period and black 
hole mass of GRS 1915+105 to be 33.5 day and IA+AMq. Interestingly, we can evolve a 
system with properties very close to this by simply starting from V4641 Sgr and following 
its future evolution with conservative mass transfer {Porb oc fi"^, where /i is the reduced 
mass); allowing for 4.6Mq to be transferred from the donor to the black hole, we have 



This would give a companion mass of 1.93M0, as compared with the Greiner et al. (2001) 
mass of Mrf = 1.2 ± O.2M0. However, the mass transfer cannot be completely conservative 
because of loss by jets, etc., as evidenced by the microquasar character of this object. 
Furthermore the above is viewed as a lower limit by Greiner et al. because the 
donor is being cooled by rapid mass loss, but its mass is estimated by comparison with 
non-interacting stars. We thus believe our evolution to be reasonable. We position the 
pre-explosion period and black hole mass of GRS 1915+105 at the same point as V4641 
Sgr. Since mass transfer and widening of the orbit always occur together, the effect of this 
post-explosion evolution is to introduce a weak secondary correlation between orbital period 
and companion mass in the long-period regime, where such a correlation is not expected to 
arise from the pre-explosion evolution. 

IL Lupi: Recently over-abundances of Mg in the companion star of IL Lupi have been 
observed (Orosz 2002). In analogy with the case of the overabundances in Nova Scorpii 
(Israelian et al. 1999, Brown et al. 2000), this indicates that there was an explosion at the 
time of black hole formation in this system, in which some of the material ejected from the 
core of the helium-star progenitor to the black hole ended up on the companion. Based on 
these observations and our given efficiency Aoce = 0.2, one can start with IIMq He star 
and I.7M0 companion as a possible progenitor of IL Lupi. From the lower boundary of the 
curve with Aoce = 0.2 in Fig. ^, the period would be 0.5 days. By losing A.2Mq during the 
explosion, the binary orbit would be widened to 1.12 day. The period had to be shortened 
to 0.8 day by magnetic braking and gravitation wave radiation before the mass transfer 
started. Conservative transfer of O.23M0 from the companion to the black hole would bring 
the period from 0.8 day to the present 1.1164 day. 




(6) 
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V404 Cyg: The black hole in V404 Cyg appears to be somewhat more massive than in 
IL Lupi, so we begin with a similar mass companion, but a IOMq black hole, which would 
have a period of 0.63 day. Again, we neglect mass loss in the explosion, although a small 
correction for this might be made later. Conservative transfer of IMq from the donor to 
the black hole then brings the period to 

0.63 day fi^^^iii2^^V = 6.7 day (7) 

close to the present 6.47 day period. Here we take IIM© and O.TM© as current masses in 
V404 Cyg (Orosz 2002). The black hole in V404 Cyg seems to be somewhat more massive 
than the others in the transient sources, with the exception of that in GRS 1915+105. 
In both cases we achieve the relatively high black hole masses and periods by substantial 
accretion onto the black hole. 

GRO J1550-564: The high mass black hole in J1550-564, 1O.56M0 (Orosz et al. 
2002), is slightly less massive than the assumed black hole mass of V404 Cyg, and the 
companion is more massive than V404 Cyg with short period, 1.552 days. So, we start from 
the same initial conditions just derived for V404Cyg (Fig. [TT|), and end up with the present 
system via simple conservative mass transfer. 

Cygnus X-1: CygX-1 is usually not considered to have come from the same evolutionary 
path as the SXTs, since it is a persistent X-ray source with a much more massive donor. 
But with the discovery of objects with relatively massive donors in the SXT category, such 
as V4641 Sgr, it is worth considering the implications of our model for it. Cyg X-1 has been 
shown to have an appreciable system velocity (Kaper et al. 1999) although it may be only 
1/3 the 50 km s~^ given there, depending on the 0-star association (L. Kaper, private 
communication). The evolution of Cyg X-1 may have been similar to that of the transient 
sources, the difference being in the copious mass loss from the companion 091 star, causing 
the black hole to accrete and emit X rays continuously. If we scale to Nova Scorpii to obtain 
the initial binary separation, we find 

^ 17.8Mq ^ ^ 

^ IMMq o o V ; 

somewhat larger than the present binary separation of 40/?©. (We would obtain 38/?© if we 
scaled from the Beer & Podsiadlowski (2001) companion mass of 2.5M0 for Nova Scorpii.) 
Given uncertainties in the mass measurements, we believe it possible for Cyg X-1 to be 
accommodated in this scheme. Some sort of common envelope envelope evolution seems 
to be necessary to narrow the orbit in the evolution involving the necessarily very massive 
progenitor stars (Brown et al. 2001a). 
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4.1. Problems with the close (AML) systems 

Reconstruction of the AML binaries is more complicated, because they have lost 
angular momentum through magnetic braking and gravitational waves, so that their present 
positions as plotted in Fig. |l] are not those at pre-explosion time. As with the evolved 
companions, matter will have been accreted onto the black hole, so the black hole masses 
will be somewhat greater than just following the explosion. As noted earlier, the binaries 
with less massive companions with separation a/ at the end of common envelope evolution 
overfill their Roche Lobes. The outer part of the companion, down to the Roche Lobe 
Rl is transferred onto the He star. This mass transfer widens the orbit to i?^, possibly 
overshooting. Unless much mass is lost in the explosion when the black hole is formed, the 
Roche lobe radius is unchanged by the formation of the black hole, and corresponds to line 
III in Fig. 0. 

Brown, Lee, & Tauris (2001c) explored the evolution of ZAMS 1.25M0 stars under 
magnetic braking, gravitational waves and mass transfer to the black hole. We adapt the 
same methods to make a more detailed study of the AML. First of all we construct (Fig. 0) 
the lower limit on the companion mass for evolution in a Hubble time, giving the dashed 
line there. All binaries with companions in main sequence at the beginning of mass transfer 
must lie between the dashed line and line HI in that figure. The fact that the AMLs tend 
to lie below the dashed line implies both mass loss from the companion and accretion 
onto the black hole. Therefore, all these systems have shrunk their orbits and increased 
their black-hole mass since the formation of the black hole, by amounts that cannot be 
determined well. In Fig. |l^, however, we show where the four shortest period AMLs would 
have come from, had they lost 0.7 Mq from an initial 1.5 Mq. From our earlier discussion 
about the a/ following common envelope evolution we saw that binaries with companions 
which stayed in main sequence were favored to come from companion masses less than 2Mq, 
and from Fig. ^ we see that they would chiefiy have companion ZAMS mass greater than 
~ IMq, so that most of them would initially have periods of 0.4-0.7 days (which follows 
from the separations obtained from our eq. (|^)). In trying to understand the detailed 
evolution of the AML we begin from a binary with a 2M0 companion which just fills its 
Roche Lobe following common envelope evolution. We then follow its evolution under the 
two assumptions made in Brown, Lee, & Tauris (2001c): (1) That its time of evolution 
is always given by its initial 2Mq mass, i.e., ignoring effects of mass loss on the internal 
evolution time (dashed lines in Fig. ^ and right dashed line in Fig. P) (2) That the evolution 
of the star proceeds according to its adjusted mass (solid lines in Fig. | and left dashed 
line in Fig. 0). Since mass loss drives the companion out of thermal equilibrium, these two 
extremes bracket the outcome of a full stellar model calculation. 
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In summary, the AML systems have had the information on their post-explosion 
parameters partly erased by subsequent evolution, in a manner that we cannot undo. 
Therefore, they can only provide a crude consistency check on the mass-period relation for 
black holes in SXTs, rather than provide precise constraints. 

5. Angular momentum and its consequences for the mass and spin of the 

black hole 

It is, in general, a difficTilt and unsolved problem to calculate the angular momentum 
of a stellar core at any given time. Even if we make the usual assumption that the rotation 
is initially solid-body, and not very far away from the maximal stable rotation frequency, 
the viscous coupling between the various layers of the star as it evolves is poorly known, 
and thus it is hard to be very quantitative. The general trend, however, is that the core will 
shrink and the envelope expand. In absence of viscous coupling, every mass element retains 
its angular momentum, and hence the core spins up as the envelope spins down, setting 
up a strong gradient in rotation frequency between the core and the envelope. Viscosity 
will then act to reduce this gradient, transporting angular momentum from the core to the 
envelope, but the efficiency of this process is very uncertain (Spruit & Phinney 1998, Livio 
& Pringle 1998). 

As we noted above (section 1), in our scenario, a number of effects will increase the 
angular momentum of the core relative to a similar core of a single star: (1) during spiral-in, 
the matter somewhat inside the orbit of the secondary is spun up by tidal torques (Rasio 
& Livio 1996); (2) the removal of the envelope halts the viscous slowdown of the core by 
friction with the envelope; (3) during the post-spiral-in evolution, tidal coupling will tend 
to spin the helium star up even closer to the orbital period than was achieved by the first 
effect. This will not be a very strong effect because the duration of this phase is short, but 
it will affect the outer parts of the helium star somewhat, and this is the most important 
part (see below). 

The net result of all these effects will be that the helium star will spin fairly rapidly, 
especially its envelope. The core is not so crucial to our argument about the fraction 
of the star that can fall into the black hole, since the few solar masses in it will not be 
centrifugally supported even in quite short orbits. For the purpose of a definite calculation, 
we therefore make the following assumptions: (1) the helium star co-rotates with the orbit 
before explosion and is in solid-body rotation; (2) the mass distribution of the helium star 
with radius is given by a fully radiative zero-age helium main sequence star. This latter 
approximation is, of course, not extremely good. However, what counts is the angular 
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momentum as a function of mass, so the fact that the mass distribution has changed from 
hehum ZAMS to explosion would be entirely inconsequential if no redistribution of angular 
momentum had taken place in the interim. As we saw above, any redistribution of angular 
momentum would take the form of angular momentum transport toward the outer layers. 
This means that relative to our ideal calculations below, a better calculation would find 
more angular momentum in the outer layers, and therefore somewhat smaller black hole 
masses than the ones we calculate. 

We now investigate how much mass will be prevented from falling into the black hole 
by the angular momentum of the He star, under the above assumptions of solid-body 
rotation with a period equal to that of the binary. If we assume that angular momentum is 
conserved during the collapse, we can get the cylindrical radius Rc within which matter is 
not centrifugally prevented from falling into the black hole: 

Rln = m^ (9) 

where /(a) is the dimensionless specific angular momentum of the marginally bound orbit 
for a given Kerr parameter a, and Mc is the total mass inside the cylinder of radius Rc. 
The Kerr parameter becomes 

where Ic is the total moment of inertia inside the cylinder of radius Rc, Ic = k'^McR^. Mc 
gives an estimate of the final black hole mass. Combining these relations with a profile 
of angular momentum and mass versus radius using the assumptions listed above, we can 
calculate the expected black hole mass and Kerr parameter as a function of SXT period 
before explosion. 



In Fig. [ly we show the predicted relation between orbital period and black-hole mass 
for different helium star masses in our model. We compare these with the present properties 
of all SXTs for which the required parameters are known. The properties are consistent 
with the theoretical relations, but do not confirm it very strongly due to the evolutionary 
changes discussed in section ^. Specifically, the AML systems lie above and to the left of the 
curves, because their orbits shrunk and their black holes accreted mass since the formation 
of the black hole. However, as we saw in section |4.1|, plausible amounts of conservative mass 



transfer since the explosion would place the systems among the theoretical post-explosion 
curves (indicated by the open squares and arrows). 

To test the theory more strongly, we show in Fig. |Tl] only those systems for which the 
pre-explosion properties could be reconstructed (section |). We compare the observed points 
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with ideal polytropic helium stars of TMq and 11 Mq, and with a full model calculation 
obtained from Woosley (2001). By coincidence, the curves converge near the region of the 
shortest-period observed systems, so that the uncertainty in helium star mass is not of great 
importance to the outcome. A helium star mass in the lower end of the range (7-9 M©) 
may be somewhat preferred for these systems. For periods above 1 day, angular-momentum 
support is not important, and the mass of the final black hole will be very close to that 
of the helium star, and thus varies somewhat from system to system. As we can see, the 
reconstructed pre-explosion properties lie much closer to the theoretical predictions. 

As a corollary, we find that systems with very large velocities, like Nova Sco, will be 
rare: at the shortest pre-explosion orbits, where much mass is ejected, the companion mass 
tends to be small. Then the center of mass of the binary is close to that of the helium 
star, which strongly limits the systemic velocity induced by the mass loss. On the other 
hand, for the widest systems, where the companion tends to be massive enough to allow a 
significant systemic velocity induced by mass loss, the mass loss itself becomes too small to 
induce much of a systemic velocity. 

An important result for our proposed relation between SXTs and hypernovae and 
gamma-ray bursts is shown in Fig. This figure shows the expected Kerr parameter of 
the black hole formed in our model. We see that for the short-period systems, this Kerr 
parameter is very large, 0.7-0.9. This means that we are justified in adding only the mass 
that immediately falls in to the black hole, because as soon as the rapidly rotating black 
hole is formed, it will drive a very large energy flux in the manner described by Brown et 
al.(2000). This both causes a GRB and expels the leftover stellar envelope. The systems 
with longer orbital periods do not give rise to black holes with large Kerr parameters, and 
thus are presumably not the sites of GRBs. 



6. Conclusions 

We have shown that there is an observed correlation between orbital period and 
black-hole mass in soft X-ray transients. We have modeled this correlation as resulting 
from the spin of the helium star progenitor of the black hole: if the pre-explosion orbit 
has a short period, the helium star spins rapidly. This means that some part of its outer 
envelope is centrifugally prevented from falling into the black hole that forms at the core. 
This material is then expelled swiftly, leading to a black hole mass much less than the 
helium star mass. As the orbital period is lengthened, the centrifugal support wanes, 
leading to a more massive black hole. The reason for swift expulsion of material held up 
by a centrifugal barrier is the fact that black holes formed in our scenario naturally have 
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high Kerr parameters (Fig. 12). This imphes that they input very high energy fluxes into 
their surrounding medium via the Blandford-Znajek mechanism, and thus power both a 
gamma-ray burst and the expulsion of the material that does not immediately fall in. 

However, because the correlation is induced between the orbital period before explosion 
and the black-hole mass, its manifestation in the observed correlation between BH mass 
and present orbital period is weakened due to post-explosion evolution of the binaries. We 
therefore considered the evolution in some detail, and for a subset of the systems were able 
to reconstruct the pre-explosion orbital periods. The correlation between pre-explosion 
period and black hole mass (Fig. |Tl|) is in much better agreement with our model than 
the original one between present period and black hole mass (Fig. |I|). We developed a 
quantitative model for the relation between period and mass, and showed that it flts the 
subset of reconstructible SXT orbits. 

Nova Scorpii stands out as the most extreme case of mass loss, nearly half of the total 
system mass, and, therefore, a great widening in the orbit which gets its period well beyond 
the gap between shrinking and expanding orbits. From Fig. ^ we see that its black hole 
mass is far below the polytropic line for its Mhg = IIM© progenitor. We believe that in 
the case of this binary a short central engine time of several seconds was able to furnish 
angular momentum and energy to the disk quickly enough to stop the infall of some of the 
interior matter not initially supported by centrifugal force; i.e., the angular momentum was 
provided in less than a dynamical time. In other words, the Blandford-Znajek mechanism 
that drives the GRB not only expelled the matter initially supported for a viscous time by 
angular momentum, but actually stopped the infall within a dynamical time. 

Since we can also compute the Kerr parameters of the black holes formed via our 
model, we find that the short-period systems should have formed black holes with Kerr 
parameters in the range 0.7-0.9. This makes them prime candidates for hypernovae and 
gamma-ray bursts, and thus provides further support for our earlier study in which we 
posited that SXTs with black-hole primaries are the descendants of gamma-ray bursts. We 
can now also refine this statement: SXTs with short orbital periods before the formation of 
the black hole have given rise to a GRB in the past. 
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Table 1: Parameters of black hole binaries in our Galaxy with measured mass functions. 
Binaries are hsted in order of increasing orbital period. All systems except Cyg X-1 (steady 
X-ray source) arc soft X-ray transients. XN indicates X-ray Nova. Earlier observations 
(Greene et al. 2001) gave the black hole mass in Nova Scorpii as 6.3 ± 0.5Mq New analyses 
of the light curve by Beer & Podsiadlowski (2001) give a somewhat smaller mass 5.4±O.3M0 
and 1.45 ± O.BSAf© for the companion. References: McClintock et al. 2001, Wagner et 
al. 2001, '^^ Bailyn et al. 1998, Filippenko et al. 1999, Gelino et al. 2001a, Harlaftis 
et al. 1996, ^) Fihppenko & Ghornock 2001. Gelino et al. 2001b, Orosz et al. 1998, 
Orosz 2002, Orosz et al. 2002, Beer & Podsiadlowski 2001, Orosz et al. 2001, 
Herrero et al. 1995, Shahbaz et al. 1994, Shahbaz et al. 1996, Greiner et al. 2001. 
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Parameter Nova Scorpii V4641 Sgr 



Orbital period (days) 
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Black hole mass (Mq) 
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Orbital separation (Rq) 


15.2 


01 QQ+l-25 
Z1.00_i 02 


Companion radius (Rq) 
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Table 2: Parameters for Nova Scorpii (Beer & Podsiadlowski 2001) and V4641 Sgr (Orosz 
et al. 2001). 
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Fig. 1. — Black hole mass as a function of present orbital period of 14 SXTs. Note that 
the orbital period is on a logarithmic scale. SXTs with subgiant or giant companions are 
indicated by big open circles (denoted as "Nu" for nuclear evolution). Filled squares indicate 
SXTs with main-sequence companions (denoted as "AML" for angular momentum loss). The 
vertical dotted line is drawn to indicate the possible existence of different classes according 
to evolutionary path of the binary, as discussed in section ^. 4U 1543—47 is marked with 
both symbols, since we believe it to be right on the borderline between main-sequence and 
evolved; for the purpose of modehng, it can be treated as evolved. 
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Fig. 2. — Companion mass as a function of present orbital period of 13 SXTs. XTE 1859+226 
is not included because the companion mass is not well determined (Filippenko & Chornock 
2001). Symbols of SXTs are the same as in Fig. ^ Line III indicates the orbital period for 
which a companion of that mass fills its Roche lobe on the ZAMS. No system can exist above 
and to the left of this line for a significant duration. Lines I and II are the upper period limit 
for systems that can come into contact while the donor is on the main sequence. For high 
masses (line II) this limit is set by the period where the evolution time of the companion is 
too short to allow the orbit to shrink significantly before it leaves the main sequence. For low 
masses (line I), where the donor never evolves off the main sequence within a Hubble time, 
the limit is set by period for which the shrinking time scale of the orbit equals the Hubble 
time. The dot-dashed line indicates the point where a system that starts its life on lines I/II 
comes into Roche contact. For very low masses, this equals line HI, because the donor never 
moves significantly away from its ZAMS radius, whereas for very high masses it equals line 
II, because the orbit cannot shrink before the companion evolves off the main sequence. At 
intermediate masses, the companion expands somewhat while the orbit shrinks, and fills its 
Roche lobe at a larger period than line HI. Systems that become SXTs with main-sequence 
donors within a Hubble time must start between line HI and line I/II. At the start of mass 
transfer, they must lie in the narrow strip between line HI and the dot-dashed line. 
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Fig. 3. — Radial expansion at the end of the giant branch. The data taken from Schaller et 
al. (1992). Radial expansions are 842i?o, 831i?0, and 750i?o for 20Mq, 25Mq, and 40Mq, 
respectively. The radial expansion shows an almost linear dependence on the ZAMS mass. 
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Fig. 4. — Possible initial orbital separations for Case C mass transfer with I.OIM© (Nova 
Sco) and 6.53Mq (V4641 Sgr) main- sequence companions. For Nova Sco, the estimated 
companion mass I.QIM© during the common envelope evolution is used as discussed in 
section ^. The x-axis is the ZAMS mass of the black hole progenitor. 
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Fig. 5. — Post-explosion orbital periods of V4641 Sgr after BH formation for various common 
envelope efficiencies. The x-axis is the ZAMS masses of black hole progenitor. The width of 
each band is determined by the initial band of possible Case C mass transfer given in Fig. ^ 
Actce ~ 0.2 is consistent with current observations. 
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Fig. 6. — An enlargement of the central region of Fig. H, indicating the present location of 
IL Lupi. If we adopt the same value Xa^e = 0.2 for the common-envelope parameter, we see 
that the reconstructed post- spiral- in period places the system right on the main-sequence 
line, implying a fairly tight constraint on the initial parameters of this system. (We took the 
ZAMS mass of the black hole progenitor to be SOM©, which corresponds to Mp ~ 25Mq in 
the beginning of case C mass transfer, and M^e = HMq.) 
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Fig. 7. — Postexplosion orbital periods of Nova Sco after BH formation for various 
common envelope efficiencies. For our calculation, we use a post-explosion black hole mass 
Mbh ~ 4.94M0 and a companion mass ~ I.QIMq. The x-axis is the ZAMS mass of 
black hole progenitor. The width of each band is determined by the initial band of possible 
Case C mass transfer separations given in Fig. |[ A post-explosion period of about 1.5 day 
is consistent with Xoce ~ 0.2. By conservative mass transfer of O.46M0 from the companion 
star to the black hole, the orbital period can be evolved to the currently observed orbital 
period. 
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Fig. 8. — Evolution of a binary with 7Mq black hole and 2Mq companion for the initial 
period of 0.54 day. The solid line marks the evolution in case the companion star adjust itself 
as it loses mass; the dashed line traces the evolution in case the mass loss does not affect 
the internal time scale of the companion star, so that it follows the same time evolution as 
an undisturbed 2Mq star. 
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Fig. 9. — Evolutionary tracks of a binary with 9Mq black hole and 2Mq companion for the 
initial period of 0.54 day. The two evolution possibilities are as in Fig. ^. Left (right) dashed 
line corresponds to the solid (dashed) lines in Fig. |[ 
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Fig. 10. — Present orbital period vs. black hole masses of SXTs. The deviations from the 
theoretical curves are substantial due to the post-explosion evolution of the binaries. The 
arrows on the AML systems point to an approximate post-explosion location if the donor 
mass was initially I.SMq, and 0.7Mq has now been transferred to the black hole. The solid 
lines indicate the possible ranges of black hole masses with polytropic index n = 3 (radiative), 
for given pre-explosion spin periods which are assumed the same as the pre-explosion orbital 
period. Here we used Rhs = O.22(Mj:/e/^0)°'^-R0- For comparison, the results with a 
"scaled" He core {9.15Mq) of Woosley's 25Mq star at the beginning of ^^C burning with 
Tc = 5 X 10^ K, appropriate for Case C mass transfer, are plotted as a dashed line (Woosley 
2001). In this plot, we scaled the radius of Woosley's core, Rwoosiey ~ 3 x 10 cm, by a factor 
2. As can be seen, the AML systems can plausible originate from systems within the curves, 
and thus are consistent with our theory. However, since they could have originated anywhere 
between the open square and their current location, they do not strongly test the theory. 
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Fig. 11. — Reconstructed pre-explosion orbital period vs. black hole masses of SXTs with 
evolved companions. The reconstructed pre-explosion orbital periods and black hole masses 
are marked by filled circles, and the current locations of binaries with evolved companions 
are marked by open circles. 
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Fig. 12. — The Kerr parameter of the black hole resulting from the collapse of a helium star 
synchronous with the orbit, as a function of orbital period. The conditions are the same as in 
Fig. |lO|, as is the meaning of the three curves. Woosley's helium core is of mass 9.15Mq from 
a ZAMS 25Mq star (Woosley 2001). Note that the result depends very little on the mass 
of the helium star, or on whether we use a simple polytrope or a more sophisticated model. 
The plot illustrates that rapidly rotating black holes needed for powering GRBs originate 
only from originally short-period SXTs. 



